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Abstract 

In this paper, we study the relationship between a diffusive model and a non-diffusive 
model which are both derived from the well-known Keller-Segel model, as a coefficient of 
diffusion e goes to zero. First, we establish the global well-posedness of classical solutions 
to the Cauchy problem for the diffusive model with smooth initial data which is of small 
L 2 norm, together with some a priori estimates uniform for t and e. Then we investigate 
the zero-diffusion limit, and get the global well-posedness of classical solutions to the 
Cauchy problem for the non-diffusive model. Finally, we derive the convergence rate of 
the diffusive model toward the non-diffusive model. It is shown that the convergence rate 
in L°° norm is of the order O (^ei^j . It should be noted that the initial data is small in 

L 2 -norm but can be of large oscillations with constant state at far field. As a byproduct, 
we improve the corresponding result on the well-posedness of the non-difussive model 
which requires small oscillations. 

Key Words: Conservation laws, chemotaxis, large amplitude solution, convergence rate, 
zero diffusion limit. 
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1 Introduction 

In this paper, we are interested in a system of conservation laws arising in chemotaxis 

pf - V • (p £ <f) = DAp £ , 

(1.1) 

qf + V (e|q e | 2 — p s ) = eAq e , 
with initial data 

£ ,q e ) 0,0) = (Po(x),m(x)) ->• 0oo,0) as |z|->-oo. (1.2) 

The chemotaxis model was preoposed by Keller and Segel in [10] to describe the traveling 
band behavior of bacteria due to the chemotactic response observed in experiments [H[2]. 
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The following Keller-Segel model has been extensively studied 



( u t = V • (DVu - xuV4>{c)) , 

(1.3) 

t Tc t = eAc + g(u,c), 

where u(x,t) and c(x,t) denote the cell density and the chemical concentration, respectively. 
D > is the diffusion rate of cells (bacteria) and e > is the diffusion rate of chemical 
substance, r > is a relaxation time scale and x > corresponds to attractive chemotaxis. 
Here g(u, c) is a kinetic function and 0(c) denoting a chemotactical sensitivity function. With 
different choices of g(u,c) and 4>(c), many results have been established in the literatures, 

cf. pa EI Eg. 

As in [HI HH H7j, if we consider the model (ll.3p with r = 1, 4>(c) = lnc, g(u, c) = —auc, 
the resulting model reads 



r ut = V • (DVu - \uV In c) , 
[ ct = eAc — auc. 

The model (jl.ip is derived from (jl.4p through the Hopf-Cole transformation 



(1.4) 



and scalings 



Vc 

q e = = - Vine, p £ = u (1.5) 

c 

t = at, x = x x j — , q = q-\/— , = — , e = — , (1-6) 

V X V « X X 

where tilde has been dropped. When the diffusion of chemical substance is so small that it 
is negligible, i.e, e — > + , then the model (jl.4j) becomes 

f ut = V ■ (DV-u - x^V In c) , 

(1.7) 

{ c t = —auc. 

A version of system (jl.7p was proposed by Othmer and Stevens in [19] to describe the chemo- 
tactic movement of particles where the chemicals are non-diffusible. The models developed 
in [19] have been studied in depth by Levine and Sleeman in [11]. They gave some heuristic 
understanding of some of these phenomena and investigated the properties of solutions of a 
system of chemotaxis equation arising in the theory of reinforced random walks. Y. Yang, 
H. Chen and W.A. Liu in [25] studied the global existence and blow-up in a finite-time of 
solutions for the case considered in [TT], respectively. For the other results on (|1.7p . please 
refer to \12\ [26] and references therein. 

Similar to the derivation of system (jl.ip , the system (II. 7p can be converted into a system 
of conservation laws as follows: 

f Pt - V • (pq) = DAp, 

(1.8) 

{ dt-Vp = 0, 

with initial data 

(p,q) (x,Q) = (po(x),q (%)) (PdojO) as |x| ->■ oo. (1.9) 

System (II. 8p has been studied by several authors. For one dimension, the global well- 
posedness of smooth solution was obtained in [27[ [6] with small initial data and large initial 
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data, respectively. For high dimensions, the global well-posedness of smooth solution to 
(|1.8p was investigated in |13l E] for Cauchy problem and initial-boundary value problem, 
respectively, where the initial data is required to be small at least in H 2 norm. For other 
related results, such as nonlinear stability of waves in one dimension and so on, please refer 
to [13 QH E3 [281 129] and references therein. 

Formally, the system (jl.ip becomes (ll.8h when we take e = 0. In fact, the investigation of 
the problem of the zero viscosity limit is one of the challenging topics in fluid dynamics and 
has been much more extensively investigated for many other models, cf. [H 02 l2~Tj l22j I24j . 
However, to our knowledge, there are few results on the system (jl.ip in this direction, cf. 
|20j . Our aim here is to prove accurately that the solutions of (jl.8p converge to the solutions 
of (jl.ip as the chemical diffusion e goes to zero. 

To do this, we first establish the global well-posedness of classical solutions to the Cauchy 
problem for the diffusive model (jl.ip with smooth initial data which is of small L 2 norm. 
Some a priori estimates independent of t and e are also obtained. Then, based on these 
estimates, we get the global existence of classical solutions to the Cauchy problem for the 
non-diffusive model (jl.8p after passing to the limits e — > 0. Finally, we derive the convergence 
rate of the diffusive model toward the non-diffusive model. 

Before stating the main results, we explain some notations. 

Notations: L p = L P (Q) (1 < p < oo) denotes usual Lebesgue space with the norm 

\\f\\ L p(n)= (jjf(x)\ p dxy , l<p<oo, 

/ = Irs 

H (Q) (I > 0) denotes the usual Zth-order Sobolev space with the norm 

ii/ii* = (eii^/ii 2 

where Q = i? 3 , and || • || = || • ||o = II ■ IIl 2 - 

The main results in this paper can be stated as follows: 

Theorem 1.1 For given number M > (not necessarily small), assume that initial data 
Oo,qo) satisfy 

Vxq = 0, ||Vpo|li2 + ||V-q ||| 2 <M, (p - Poo, qo) G H 3 , Poo > 0. (1.10) 

Then there exists a positive constant eq depending on M andpoo such that the Cauchy problem 
exists a unique global solution in R 3 x (0, oo), which satisfies 

' (P £ -Poo,q £ ) G L°°([0,oo),F 3 ), 

< Vp £ G L 2 ([0,oo),F 3 ), Vq £ G L 2 ([0, oo), F 2 ), (1.11) 
k eVq £ G L 2 ([0,oo),F 3 ) 
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and 

l|q £ WII^ + \\ P £ (t)- Poo \\ 2 H3 



1.12) 



+ f (\\Vf(s)f HS +e\\Vcf( S )f H3 )ds+ [ t \\Vcf(s)\\ 2 H2 ds<C, 
Jo Jo 

where C is a positive constant independent of e and t, provided that 

||qoll£ 2 + ||po - Poo\\ 2 L 2 < e . (1.13) 

The last theorem is concerned with the convergence rate as well as the global well- 
posedness of (Oft - lfO!) . 

Theorem 1.2 Suppose that (po, qo) satisfies the assumptions in Theorem ! 1 . 11 then the Cauchy 
problem \1.8\ )- [T79\) exists a unique global solution in R s x (0, oo), which satisfies 

f (p-Poo,q)G^ oo ([0,oo), J ff 3 ), 

(1.14) 

[ Vp£ L 2 ([0,oo),F 3 ), VqeL 2 ([0,co),ff 2 ). 

Furthermore, 

W-v) ®\\m + ll(q £ - q) < (i-is) 

In particular, 

IK/ -P) + ||(q £ - q) (*)|U < Cei (1.16) 

Here C is a positive constant independent of e and t. 



Remark 1.3 Notice that for the global existence of the solutions to Cauchy problem il.ty) - 
il.9\) , we only assume that the L 2 norm of initial data is small. The initial data can be of 
large oscillations with constant state at far field. This is an improvement of ]13\j where the 
initial data is required to be small in H s (s > | + 1) norm which implies the oscillations are 
small. 

Remark 1.4 The power of e in jil.l5\) could be improved to 1, which needs a slightly modi- 
fication of the proof of Theorem with more regular initial data. But in this case, it seems 
that the coefficient C might depend on t. 

The proofs of Theorems 11.11 and 1 1 . 2 1 are based on the classical energy method. The key point 
for the proof of Theorem 11.11 is to obtain some a priori estimates independent of e in which 
the L 2 — bound of V • q £ plays a crucial role. For the proof of Theorem II .21 some estimates of 

the order O are required, which needs some delicate analysis. 

The rest of the paper is organized as follows. In Section 2, we study the global unique 
solvability on the Cauchy problem (|1.1|) - (|1.2|) . In Section 3, the zero-diffusion limit as well 
as the global well-posedness of the solutions to (jl.8p - (jl.9p is considered. We show that the 
convergence rate in L°°-norm is of the order 0(e^), when diffusion parameter e — > + . 

Throughout this paper, we denote a generic positive constant by C which is independent 
of e and t. 
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2 Proof of Theorem 11.11 



In this section, we are concerned with the global existence of large-oscillations solutions to 
the Cauchy problem (jl.ip - (jl.2p when the initial data is sufficiently close to a constant in 
L 2 -norm. The global existence follows from a local existence theorem and some a priori 
estimates globally in time. 

The local existence of the solutions could be done by using some arguments similar to 
[13j . We shall get some a priori estimates globally in time which are also uniform for e. 

More precisely, for any given T > and e > 0, suppose (p £ (x, t), q £ (x, t)) is a smooth 
solution to the Cauchy problem (|l.ip - (|1.2[) with regularities as in Theorem II. 1\ we get the 
following key proposition. 

Proposition 2.1 For given number M > (not necessarily small), assume that initial data 
(Po,qo) satisfy 

Vxq = 0, HVPollia + llV-qollln < M, (po " Poo, q ) G H 3 , Poo >0. 

Then there exists some positive constant Eq depending on M andpoo such that if(p £ (x, t), q e (x, t)) 
is a smooth solution to the Cauchy problem of 111. 1]) - [TJ^) in M. 3 x (0, T], satisfying 

||q £ || 2 2 <2e , ||V-q e || 2 2 < 2M, (2.1) 

then 

||q £ |l! 2 <^ , ||V-q £ ||| 2 <^M, (2.2) 

provided that 

llqollia + ||po -PooWl* < ^0- 

Without loss of generality, we suppose that D = l,Poo = 1- Letting p = p — 1, we obtain 
that 

f Pt ~ ^ ' (p e q e ) — V • q e = Ap £ , 

{ „ ( 2 - 3 ) 

I q| + V (e|q £ | 2 - p e ) = eAq e , 

with initial data 

(F,q e )(x,0) = (po-l,q ) -> (0,0) as \x\^ao. (2.4) 



Proof of Proposition I2.lt 

The proof of Proposition 12.11 consists of the following Lemmas 12.21 and 12.31 
Lemma 2.2 [1? estimate] Under the conditions of Theorem \l. 11 it holds that 



|~£||2 

\p II 



+ llq £ f + T WWW 2 ds + e f ||Vq e || 2 ds < ^£ , (2-5) 

JO JO 2 



provided that eq is small enough. 

Proof. Multiplying the first equation in (|2.3p by 2p e and the second by 2q e , summing up 
them and then integrating over M 3 x [0, t] , one gets after integration by parts that 



+ \\q £ \\ 2 + 2[ \\Vp^\\ 2 ds + 2e [ ||Vq £ f ds 
Jo Jo 

+ 1 1 qg 1 1 2 + 2 J J p^cf ■ (Vp^)dxds + 4e J J (q £ ) T • (Vq £ • cf)dxds. 



(2.6) 
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Next, we shall estimate the last two terms in the right-hand side. By Cauchy inequality, 
Holder inequality, Sobolev inequality and Gagliardo-Nirenberg inequality, we obtain 

2/ /Vq E • (W)dxds < [ \\Wf ds + C f Wftffds 
Jo J Jo Jo 

< fwWW 2 ds + C f\\n%W\^ds 
Jo Jo 

<[ WWW 2 ds + C f \\VF\\ 2 \\Vtf\\\\tf\\d8 (2.7) 
Jo Jo 



and 



4e f* I (q E ) T • (Vq E • cf)dxds < 4e f ||q e || L3 ||Vq E || L2 ||q ( 
Jo J Jo 



'" L e ds 



<Ce I ||Vq E ||t||q e ||f ||Vq £ || 2 ds. (2.8) 
Jo 

Since Aq e = V(V • q £ ) - V x (V x q e ), we obtain that 

Vq e = -V(-A)- 1 V(V • q e ) + V(-A)- X V x (V x q E ). 
The standard L 2 estimate shows that 

||Vq £ || i2 <C(||V-q £ || L2 + ||Vxq £ || L2 ). (2.9) 
Moreover, by taking the curl for (12.3p ^>. one has 

^(Vxq £ ) = eA(Vxq e ). (2.10) 

Initial data is given as 

(V x q £ ) | t=0 = V x q = 0. (2.11) 
By solving the initial value problem (|2.1UI ) and (|2.11l ) , one has 

V x q e = 0, (2.12) 

which implies 

Aq e = V(V-q e ), ||Vq £ || i2 < C||V-q £ || i2 . (2.13) 
The combination of flU]), (pT7|) . and (f2T3l) yields 

2/ /Vq £ ■ (W)dxds < [ WWW 2 ds + C^/mTq f ||VF|| 2 ^, (2.14) 
Jo J Jo Jo 

and 

4e [ /'(q £ ) T -(Vq £ -q £ )(ixds<Ce(Meo)3 I ||Vq £ || 2 ds. (2.15) 
Jo J Jo 

Substituting (|2.14p and (|2.15p into (|2.6p and setting £0 < one may arrive at (|2,5p . 

where (|1.13p has been used. This completes the proof of Lemma 12.21 □ 
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Lemma 2.3 [First-order energy estimate] Under the conditions of Theorem Xl.lX it holds that 
WWf + l|V • q £ f + ^ ||V • q £ || 2 + J* f ds + e || Aq £ f) ds < ™, (2.16) 

provided that eq is small enough. 

Proof. Notice that 

V • qf = Ap 6 - eA|q £ | 2 + eA(V • q £ ) 

= Pi ~ V ■ (p £ q £ ) - V • q £ - eA|q £ | 2 + eV • (Aq £ ), (2.17) 

where we have used f|2.13|> . 

Multiplying (|2.17p by 2V • q £ and integrating the resulting equation over M 3 , one obtains 
after integration by parts that 

4||V-q £ || 2 + 2||V-q £ || 2 + 2e||Aq £ || 2 
at 

= 2 J gV • cfdx - 2 J V • (gq £ ) V • tfdx + 2e J V|q £ | 2 (Aq £ )dx. (2.18) 

Next, multiplying the first equation in (j2.3|) by 2p £ , integrating the resulting equality over 
R 3 and using integration by parts, one has 

|||Vif|| 2 + 2||i5f|| 2 = 2 |i5fV-q £ cte + 2 J V ■ (jfcf )jf t dx. (2.19) 

The combination of f[27T8|> with (|2~T9j) yields 



(||V • q £ f + \\W\fj + 2 ||V • q £ f + 2 ||pf || 2 + 2e || Aq £ || 2 

= 4 y pfV • q £ dx + 2 f^- (j?<fWtdx 

-2 J V • (pV) V • q £ dx + 2 ^ y V|q £ | 2 ■ (Aq £ )dx 
4 

= ^J,. (2.20) 

i=l 

For Ji, using (|2.17p . and integration by parts, we have 
Ji = 4-^ y V • tfjFdx - 4 y V • qfgiix 

= 4 HVpl 2 + 4-^ y V • q e pFda; - 4e y V|q £ | 2 • (VpF)dx - Ae J Aq £ • (Vp^)dx 

3 

= 4||Vp E || 2 + ^ J|. (2.21) 



i=i 



By Cauchy inequality, Holder inequality, Sobolev inequality, Gagliardo-Nirenberg inequality 
and (|2.ip . we estimate J^-Jf as follows: 

J\ = -8e y q £ • (Vq £ ) • (Vp^)dx 
<C6||q £ || L 3||Vq £ || L6 ||Vp £ || L 2 

< Ce(Me )^ ||Aq £ || ||Vp^|| 

< Ce (Me )^ || Aq £ || 2 + C WWW 2 (2-22) 



and 



Jl< i E ||Aq 6 f + CUV? 6 !! 2 . (2.23) 
On the other hand, Cauchy inequality gives 

J 2 + J 3 < \ ||V • q £ || 2 + \ W t f + 16 ||V • (?f q £ )|| 2 • (2.24) 
Similar to (|2.22p . it is immediate to obtain 

J 4 < Ce{Msq)\ ||Aq £ || 2 . (2.25) 
Finally, we estimate the last term on the right hand side of (|2.24|) 

6 

II V • (pFq e )|| 2 < C ||VpF • q £ || 2 + C ||if V • q £ || 2 = £ Jj. (2.26) 

i=5 

For J5 and J6, using Sobolev inequality, Holder inequality, Gagliardo-Nirenberg inequality, 
Young inequality, (|2.ip and (|2. 13|) . we obtain 

J5<C||VF|||6||q E ||i3 

<C||V 2 F|| 2 ||Vq e || ||q £ || < C^/MT \\V 2 ^\\ 2 (2.27) 

and 

J6<c(\m\l* + \\VF\\l*) liv-q £ ll 2 

< C (||pF||^ HV^Hi + ||Vp £ ||5 HV^p) ||V • q £ || 2 
<C||Vif|| 2 ||V-q £ || 2 + Ce ||V-q £ || 2 

+C||Vp £ || 2 ||V-q £ || 8 + ei ||V 2 pP|| 2 

< CM(1 + Af 3 ) IIV^H 2 + e x IIV 2 ^)! 2 + Ce ||V • q £ || 2 . (2.28) 
Substituting (p^7D - (l2T28l ) into (i!T2"6]) shows that 

II V • (Fq e )H 2 < (c^/mFv + ei ) IIV 2 ^!! 2 + CM(1 + Af 3 ) llV^f + Ce ||V • q £ || 2 . (2.29) 

Applying the standard L 2 -estimate for (j2.3j) i . one has 

||v 2 if || 2 < C (||V • q £ || 2 + ||g|| 2 + ||V • (?f q £ )|| 2 ) , (2.30) 

which together with (12. 29ft gives 



II V • (F q £ )|| 2 < [C^/mTo + CeiJ (||V • q £ || 2 + m\\ 2 + IIV • (jf q £ )| 

+CM(1 + Af 3 ) || VpF|| 2 + Ce ||V • q £ || 2 . (2.31) 
Taking ex suitably small, one can deduce that 



1 ~e ri2 



||V • (if q £ ) || 2 < (C^/M^+Cex) (||V-q 6 |f + NW | 

+CM(1 + Af 3 ) ||Vpl 2 + Ce || V • q £ || 2 , (2.32) 



ds 



, ^ (1 - 2Cei) 2 
provided £o < 4C2M • 

Substituting (|2.2ip - (|2.25j) and f|2.32[) into (|2.20[) and integrating the resulting inequality 
over [0, t), one gets after using (jl.lOp . (|1.13|) and (|2.5p that 

l|V • q e f + ||VFf + 2 jT ||V • q e || 2 ds + 2^ (||pf f + e ds 
<M+ (cy/M^ + Ce x + Ce + ^ J q (ll V " tff + Mil*) ds 

+Ce (M + M 4 + l) + ^C(Meo)^ + ~) e ^ || Aq e || 2 

+4 y V • cfjfdx - 4 y V • qoPocte, 
which together with Cauchy inequality and (|2.ip deduces 

||V • q e || 2 + llVFf + 2 £ ||V • q e || 2 ds + 2 £ (||g|| 2 + e || Aq £ || z ] ds 
< ™± + (^Cx/M^ + Cei + Ce + ^ f (|| V • q £ || 2 + ||g 

+Ce (M + M 4 + l) + ^C(Meo)^ +0 e ||Aq e || 2 ds. (2.33) 

Next, choosing s\ suitably small and taking 

. f 1 M 1 - 2C £l 1 

eo _ mm | , + M ± + 1 y 2C( ^ + 1} ' J ' 

one can get (|2.16p . This completes the proof of Lemma 12.31 □ 
The proof of Proposition 12.11 is complete. To prove Theorem 11.11 we need to get some 
high order estimates. Before beginning, we give the following corollary. 

Corollary 2.4 Under the conditions of Theorem \l.l\ it holds that 

l\\\V 2 m 2 + \¥\\l™)ds<C, (2.34) 
J o 

provided that £q is small enough. 

Proof. It follows from (HOU!) . (12311 and (I2T61) that 

r||V 2 p E || 2 ds<C'fr||V-q e || 2 ds+ f \\p £ t \\ 2 ds + f \\V ■ (^q E ) || 2 ds) 
Jo \Jo Jo Jo J 

< C. (2.35) 

By Sobolev's embedding theorem, combining (|2.5p and (I2.35p . we get (I2.34j) . □ 

Lemma 2.5 [Second- order energy estimate] Under the conditions of Theorem it holds 
that 

||V 2 p £ || 2 + ||V(V • q e )|| 2 + f ||V(V • q £ )|| 2 ds 

Jo 

+ f (ll Wt II 2 + e || V • (Aq e )|| 2 ) ds < C, (2.36) 
provided that £q is small enough. 
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Proof. Differentiating (|2.17p yields 

V(V • qf) = V(Ap e ) - eV(A|q £ | 2 ) + eV(A(V • q £ )) 

= Vpf - V (V • CpV)) - V(V • q £ ) - eV(A|q £ | 2 ) + eV(V • (Aq £ )). (2.37) 

Multiplying ()2.37|) by 2V(V • q £ ) and integrating the resulting equation over R 3 , one obtains 
after integration by parts that 

| ||V(V • q £ )f + 2 ||V(V • q £ )|| 2 + 2e ||V • (Aq £ )|| 2 

= 2 J V(V • q £ ) • {W t )dx - 2 J V (V ■ (pFq 6 )) • (V (V • q £ )) dx 

+2e J A\cf\ 2 V ■ (Acf)dx. (2.38) 



Next, applying V to (|2.3P i . multiplying it by 2Vpf, integrating the resulting equality over 
R 3 and using integration by parts, one has 

j t || Apl 2 + 2 ||Vj5f || 2 = 2 J V(V • q £ ) • (Vj5f)dz + 2 1 V(V • (pV)) • (Vg)dz. (2.39) 
Putting (f2T38D and (1239"]) together, we have 

j t (||V(V • q £ )f + || Apl 2 ) + 2 ||V(V • q £ )f + 2 ||Vpf || 2 + 2e ||V(Aq £ )|| 2 
= 4 y -V(V • q £ ) • (Wt)dx + 2 y V(V • (Fq £ )) • (V^)dx 

-2 y V(V • (Fq £ )) • (V(V • q £ )) dx + 2e y A|q £ | 2 V • (Aq £ )dx 

10 

= ^J t . (2.40) 



i=7 



For J7, using (|2.3T|) . and integration by parts, we have 

J 7 = J V(V • q £ ) • {W)dx -A J V(V • q £ ) 4 • (VfP)ds 

= 4 HApFH 2 + 4:4- ! V(V • q £ ) • CVjF)dx 
dt J 

-As y A|q £ | 2 Ap £ dx + 4e J V • {A^Ajfdx 
3 

= 4||Ap £ || 2 + ^jf (2.41) 



For J 2 , using Sobolev inequality, Holder inequality, Gagliardo-Nirenberg inequality, Young 
inequality, Proposition 12.11 and (|2.13p . we obtain 



J2 = - 8e y (|Vq £ | 2 + Aq £ • q £ ) Ajfdx 



< - || Apl 2 + Ce (||Vq £ ||i 4 + ||q £ Aq £ || 2 

< | || A^H 2 + Ce (||Vq £ ||^ ||V 2 q £ ||^ + || Aq e ||| 6 ||q £ ||i 3 

< 1 1| A^f + Ce || Aq £ || 2 + ||Vq £ || 2 + Ce ||V(Aq £ )|| 2 ||q £ || ||Vq £ 



< - || A^H 2 + Ce ||Aq £ || 2 + Ce ||V • q £ || 2 + Ce^/Me~ ||V • (Aq £ )|| 2 , (2.42) 
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where in the last inequality we have used the following fact 

||V(Aq £ )|| 2 < C (||V • (Aq £ )|| 2 + ||V x (Aq £ )|| 2 ) 

= c(||V-(Aq e )|| 2 +||A(Vxq e )|| 2 ) 

= C||V-(Aq £ )|| 2 
due to (j2.9j) and (|2.12p . Cauchy inequality gives 

J*< l -e\\V-(A<f)f + C \\AW\\ 2 (2-43) 

and 

J 8 + J9 < \ ||V(V • q £ )|| 2 + \ ||Vpf|| 2 + 16 ||V(V • (W))!! 2 . (2.44) 

Similar to J 2 , we estimate J\q as follows: 

J 10 < | ||V • (Aq £ )|| 2 + Ce^M^ ||V • (Aq £ )|| 2 

+Ce||Aq e || 2 + C7e||V-q e || 2 . (2.45) 
Finally, we estimate the last term on the right hand side of (12.44ft 
||V(V-(^q £ ))|| 2 

< C ||V 2 F • q E || 2 + C HV^V • q E || 2 + C \\W • (Vq e )|| 2 + C ||FV(V • q e )|| 2 
14 

= Yl J i- ( 2 - 46 ) 
i=ll 

By Holder inequality, Sobolev inequality, Gagliardo-Nirenberg inequality, Young inequality, 
Proposition 12.11 and (12. 13H . we estimate J11-J14 as follows: 

Jii<C||V 2 F||^||q £ |li3 

< C yV 3 ^!! 2 ||Vq £ || ||q e || < Cy/MFQ\\^f\\ 2 , (2.47) 

J12 + Jl3<C \\W\\l°o (||Vq £ || 2 + ||V • q e || 2 ) 

< ClIV 2 ^!! IIV^H ||V-q £ || 2 

< C ||V • q £ || 4 IIVVH 2 + e x \\V Z W\\ 2 (2-48) 

and 

Ji4<C||if||^||V(V. q £ )|| 2 . (2.49) 
Substituting (p^7D - (f2T4"9l ) into (12^61) shows that 
||V(V-(Fq £ ))|| 2 

< (Cy^o + e!) \\V 3 f\\ 2 + C p 2 ^\\ 2 + C W\\ 2 Laa ||V(V • q £ )|| 2 . (2.50) 
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Applying the standard i7 1 -estimate for ()2.3|) i leads to 

||V 2 if||^ < C (||V • q £ ||^ + m\W + IIV ■ (¥°f)\\ 2 m) , (2-51) 
which together with (|2.50p gives 



V(V • (Fq £ ))f < <Vm^ + Ce l )( ||V • q £ ||^ + ||pf ||^ + ||V • {ffcf)\\ 2 m 



+C ||V 2 ^|| 2 + C ||V(V • q e )|| 2 . (2.52) 

Substituting (|2.41j) - (|2.45j) and (|2.52p into (|2.40p and integrating the resulting inequality over 
[0,t), one gets after using (jl.lOp . Lemmas l2.2H2.3l and Corollary 12.41 that 

||V(V-q £ )|| 2 + ||AF|| 2 

+2 / < ||V(V-q £ )|| 2 ( fe + 2 \\X7fi\\ 2 ds + 2e f ||V • (Aq £ )|| 2 ds 
Jo Jo Jo 

< C + 1 ||V(V ■ q £ )|| 2 + £ W\\l~ ||V(V • q £ )|| 2 ds 

+e [ ||V- {Acf)\\ 2 ds + Ce^JW^Q [ ||V • (Aq e )|| 2 ds 
Jo Jo 

+ (c^JmF + Ce x + ^ £ (||V(V • q £ )|| 2 + ||Vpf|| 2 ) ds. (2.53) 

One obtains (I2.36P by using Corollary 12.41 and Gronwall's inequality. This completes the 
proof of Lemma 12.51 □ 

Corollary 2.6 Under the conditions of Theorem \l.l\ it holds that 

!\\\V 2 ¥fm + l|VF|li~) ds <C, (2.54) 
J o 

provided that Eq is small enough. 

Proof. It follows from ([23TD - (p32D . ([2TT6]) and (p36l) that 

rilV^II^ ds 
Jo 

<c(J \\V-cf\\ 2 Hl ds + \m\\mds + J ||V-(Fq £ )||^i^) <C. (2.55) 
By Sobolev's embedding theorem, we get 

/"* HV^Hioo ds<<7, 
J o 

which together with ([235]) leads to ([234]) . □ 

Lemma 2.7 [Higher-order energy estimate] Under the conditions of Theorem it holds 
that 

nt 

||V 3 p £ || 2 + ||V 2 (V-q £ )|| 2 + / ||V 2 (V-q e )|| ds 

Jo 

+ (|| v2 pf|| 2 + e ll v3 ( v -q £ )|| 2 ) ds < c, ( 2 - 56 ) 



o 



provided that sq is small enough. 
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Proof. Applying A to (|2.17p . multiplying it by A(V • q £ ), taking integrations in x and 
using integration by parts, one gets 

1 1| A(V • q e )|| 2 + 2 || A(V • q e )|| 2 + 2e || A 2 q £ || 2 

= 2 J Apf A(V • cf)dx -2 j A(V • (p E q e ))A(V • cf)dx 

+2e J V(A|q E | 2 ) • (A 2 q £ ) dx. (2.57) 

Similar to (|2.57p . one has 

I l|v(A?)l12 + 2 imf ~ 2 1 A5 < A(v • + 2 1 A(v ■ <W))A ^ (2 ' 58) 

Putting the above two equalities together, we get 
4(||A(V-q £ )|| 2 + ||V(Ap £ )|| 2 



dt 

+2 || A(V • q e )|| 2 + 2 || Apf || 2 + 2e || A 2 q £ || 2 
= 4 J Apf A(V • cf)dx + 2 y A(V • (p £ q £ ))Apfdx 

-2 / A(V • (^q £ ))A(V • cf)dx + 2e j V(A|q e | 2 ) • (A 2 q e ) dx 

18 

= J2 J i- ( 2 - 59 ) 

i=15 

In a manner similar to the estimates of J7-J10, Jib-Jid, can be bounded as follows: 

J 15 = 4— y A(V • q e )Ap e da; - 4y A(V • q e ) t AF*c 

= 4||V(Ap E )|| 2 +4-| y A(V • q e )ApFda; 

-4e y A 2 q £ - (V(Ap E ))dx-4e y V(A|q £ | 2 ) • (V(A^)) dx 
3 

= 4||V(A^)|| 2 + ^jj 5 . (2.60) 



i=l 

Cauchy inequality gives 



jf 5 <~e\\A 2 cf\\ 2 + C\\V(AF)\\ 2 . (2.61) 



For Jf 5 , we have 



Jf 5 < I ||V(Ap £ )|| 2 + (llVq^l^ || Aq £ || 2 + ||V(Aq £ )||i 6 ||ql£ 3 

< |||V(Ap £ )|| 2 + Ce||vV|| ||V 2 q E || ||Aq E || 2 + CeyM^||A 2 q e || 2 

< I \\V(A^)\\ 2 + Ce || A(V • q e )|| 2 + Ce || Aq £ || 2 + Ce^M^ || A 2 q e || 2 . (2.62) 

By Cauchy inequality, we get 

Jie + Jn < \ || A(V • q £ )|| 2 + \ ||A^|| 2 + C || A(V • (Fq £ ))f • (2-63) 
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Similar to Jf 5 , J\$ can be estimated as follows: 

Jia < |||AV||^Ce||A(V-q £ )|| 2 + C'e||Aq e || 2 + C'ev / M^||A 2 q E || 2 - (2-64) 
Next, we estimate the last term on the right hand side of (|2,63p 
l|A(V-(Fq £ ))|| 2 

< C ||V(A^) • q £ f + C ||V 2 F • (Vq e )|| 2 + C HAFV • q e f + C ||V(V • q e ) • (V^)l| 2 
+C||Vp E -(A q £ )|| 2 + C||p E A(V- q £ )|| 2 

24 

= J i- ( 2 - 65 ) 
i=19 

We estimate J19-J24 as follows: 

J 19 < C||V(AF)||| 6 ||q B ||£a < C||A 2 i?|| 2 ||Vql ||q £ || < C^M^ || A 2 ^|| 2 , (2.66) 

J20 + hx < C \\V 3 ^\\ \\V 4 ^\\ ||V • q e || 2 

< C ||V - q £ || 4 ||V(Z\p-)|| 2 + £l || A 2 ^!! 2 

<C||V(A^)|| 2 + ei ||A 2 n| 2 , (2.67) 
J22 + J23 < C IIV^H 2 ^ ||V(V • q £ )|| 2 < C \\W\\ 2 Lao (2.68) 

and 

J 2 4<C||p E || 2 00 ||A(V-q £ )|| 2 . (2.69) 
Substituting (p^6l) - (l2T69l ) into (i!T65l) shows that 

^(V-^))!! 2 ^ (Cv/M^ + ei) ||AV|| 2 + C||V(A^)|| 2 

+c \\W\\U + c \w\\U II A(V • q £ )f . (2.70) 
Applying the standard i7 2 -estimate for (|2.3P i leads to 

||A^||^ 2 < c (||v • q £ ||| 2 + Wtfm + liv • (F<f)\\h) > ( 2 - 71 ) 

which together with (|2.70p gives 

C || A(V • (i?q £ ))|| 2 < (CVM^ + Cei) (||V • q £ ||^ + llpfll 2 ,, + ||V • (pV)!!^) 

+C ||V(AF)|| 2 + C HFIlic II A(V • q £ )|| 2 + C WWWloo ■ (2.72) 

Substituting (|2.60p - (|2.65|) and (|2.72|) into (|2.59|) and integrating the resulting inequality over 
[0, t), one obtains after using (jl.lOp . Lemmas l2.2M2.31 Lemma [2.5| Corollary 12.41 and Corollary 
US] that 

||A(V-q £ )|| 2 + ||V(A^)|| 2 

+2 J* ||A(V • q e )|| 2 ds + 2 J* (||Apf || 2 + e ||A 2 q e || 2 ) ds 

<C+ l - || A(V • q £ )|| 2 + (c^/mT + Ce 1 + ~) j* (|| A(V • q £ )|| 2 + || Apf || 2 ) ds 

+ C I HpIL l|A(V-q e )|| 2 ds + e I \\A 2 q e \\ 2 da + Cey/Me^ [ ||A 2 q £ || 2 ds. (2.73) 
Jo Jo Jo 
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Using Corollary 12.41 GronwalFs inequality and (|2.13p . one can immediately get (|2.56p . This 
completes the proof of Lemma 12.71 □ 

Corollary 2.8 Under the conditions of Theorem \l.l\ it holds that 

f \\V 2 p e \\ 2 H2 ds<C, (2.74) 
J o 

provided that Sq is small enough. 

Proof. It follows from (l2~T6D . (12T36D . (12361) and (1277T]) that 

2 ^\\ 2 H2 ds<C (J ||V • cf\\ 2 H2 ds + J \\m 2 H 2 ds + j ||V • (Fq £ )llH 2 < C- 

□ 



t 

iiv^ir w2 ^ < 







Proof of Theoremll.lt 

As a consequence of ([23]) . ([235]) , (1236]) . (1236]) and flZHj) , one obtains 

l|q £ WllH3 + IIFWIlk + /" (HVfFOOIII, + £||Vq e ( S )||^) d S 

J 

+ / ||Vq £ ( S )||| 2 d S <C, 
JO 

where we have used (|2.12p and the following fact: 

llVq 6 !)^ < C(||V-q e || Hs + ||V x q £ ||^ s ), for s = l,2,3. (2.75) 

The uniqueness result can be proved by the method used in p3], we thus omit the details 
for brevity. This completes the proof of Theorem 11.11 

3 Proof of Theorem Q 

Prom (jl,12p . one can easily get a unique smooth solution (p, q) to ((1 .8[> - (|1.9[) with regularities 
as in Theorem 11.21 after passing to the limits e — > (take subsequence if necessary) . 

We will give the proof of the last part of Theorem II .21 i.e., the convergence rate as e — > 0. 
To do this, it suffices to show the following Proposition 13.11 

Proposition 3.1 Assume that the assumptions listed in Theorem li.il are satisfied. Then 
there exists a positive constant C , independent of t and e, such that 

||q e - q||#2 + \\p £ - p\\ 2 H 2 
ft 



[ ||V(q e -q)||^ + ||V(p £ -p)||^d S <C £ . (3.1) 
J o 



+e 

Proof. Denote D = 1, p = p — 1. Then (jl.8p - ([1.9p can be translated into the following 
system 

f pt - V • (pq) - V • q = Ap, 

(3.2) 

{ ca-Vp = 

with initial data 

(p,q) (x,Q) = (p (x),qo(x)) ->• (0,0) as \x\ -> oo. (3.3) 
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Setting 

ip £ = q £ - q, 9 s = p £ - p =j? -p. (3.4) 

Then we deduce from ([22])-([23]) and ([23D-([23]) that (ip £ , 9 £ ) (x, t) satisfy the following 
Cauchy problem: 

9% - V • (ij £ p + q6 £ ) - V • V e = A0 £ , 
+ V ( e (q £ ) 2 - 9 £ ) = eAifj £ + eAq 
with initial data 

(V> £ ,0 e )(x,O) = (O,O). (3.6) 

Step 1 . 

Multiplying the fist and second equations of (|3.5p by 29 s and 2^ £ respectively, integrating 
the adding result with respect x and t over M 3 x [0, i] , we have 

WW 2 + P e \\ 2 + 2 f Q ( e iiwi 2 + ll vr II 2 ) ds 

= -4e y ^ q £ • (Vq £ ) • V £ rfxds +2e J J i/) £ ■ (Aq)dxds 
-2 y y pW • ^ctecfs - 2 y y q . (V9 £ )9 £ dxds 

= J2 K i- (3-7) 

i=i 

By Cauchy inequality, Holder inequality, Sobolev inequality, Gagliardo-Nirenberg inequality 
and Theorem I Lit we obtain 



K x <2e I \\Vcf\\Us + 2e / ||q E ^ e || z ds 

(is 



<ce + G6 ri| q £ iii3iiv £ iii6 

JO 

<Ce + Cey^Me^ [ ||Vi/> £ || 2 ds. (3.8) 

JO 



Integration by parts and Cauchy inequality implies 



K 2 <e [ ||Vq|| 2 ds + e f \\V^ £ fds 
Jo Jo 



(3.9) 



and 

K 3 +K 4 <2^ (llpll^ + HqH 2 ^) (ll^f + H^ll 2 



ds 



+^ I ||Vrfds. (3.10) 



Substituting (l3T8l)- (^T0l) into §3li), we get 

Wt + WW 2 + y (e IIV^II 2 + l|V0 £ || 2 ) ds 



'o 

£ ^ - ^ / + llqllloo) (ll^f + ll^f) ds. (3.11) 
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Step 2 . 

Multiplying the fist and second equations of (|3.5p by —2A6 £ and —2Aip £ respectively, inte- 
grating the adding result with respect x and t over M 3 x [0, t] , we have 



ds 



HV-^II 2 + ||V0 E || 2 + 2 J q ( £ \\ A ^ £ f + \\ Ad£ 
4e J J cf ■ (Vq e ) • Aip £ dxds - 2e j Aip 6 ■ (Aq)dxds 
-2 J J V • (ip £ p + q6 £ ) A9 £ dxds 

7 

Y, K i- ( 3 - 12 ) 



i=5 

Cauchy inequality leads to 



and 



K b <Ce f \\cf\\l 00 \\Vc L £ \\ 2 ds + ^ f \\Aij £ \\ 2 ds (3.13) 
Jo * Jo 



K • / ||Aq £ || 2 (is + e / \\Ai/ffd8. (3.M) 



Straightforward calculations show that: 



K 7 = -2 f I V • i> £ pA6 £ dxds -2 f V e • {W)^0 £ dxds 

Jo J Jo 

-2 j j V • q(9 e A6 £ dxds -2 j J q • (W E ) A6 £ dxds 

4 

= (3-15) 



i=l 



Here if^ — iTy are estimated as follows: 



ft /-t 
K$<^l \\A9 £ \\ 2 ds + C j \\p\\ 2 LOO \\V^ £ \\ 2 ds, (3.16) 



o jo 



K%<][ \\A6 £ \\ 2 ds + C f ||Vq e ||£oo ||V e || 2 ^ (3.17) 
4 io Jo 



and 

i r l 



K% + K$<\! \\A6 £ \\ 2 ds + C f ||V-q|||oo W t ds 
4 Jo Jo 

+c ! \H\l°° livrfds. 

Jo 

Substituting ([37T3T) - (f37TH] ) into (l3~T2t we get 



||VV/f + ||V0 £ || 2 + j (e || A^ e || 2 + || A9 £ f 



o 



(3.18) 



<Ce + C I ( ||p1|^ 1i00 + ||q||^ 1)00 ) ( ||^||^ + ||0 e ||^ ) ds. (3.19) 
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Step 3. 



Differentiating (|3,5p yields 

{ V0f - V (V • (ip e p + q_6 £ )) - A^ £ = V(A6< £ ), 

V • z/f + A (e (q £ ) 2 - = eV • (AV> £ ) + eV • (Aq). 



Multiplying the fist and second equations of (pT20j) by -2V(A6> £ ) and -2V • (Aip £ 
tively, integrating the adding result with respect x and t over R 3 x [0,t], we have 



(3.20) 



respec- 



|A^ £ || + ||A6 



+ 2 



o 



£||V-(A^)|| 2 + ||V(Ar)|| 2 W'-: 



-2 y y V (V • (V/p + q# £ )) • (V (A0 e )) dxds 
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i=8 

Then, it follows from Cauchy inequality that 

2 



if 8 = J V • (A^ e ) \Vcf\ 2 dxds + te J V • (A^ £ )q e • (Aq £ )dxds 

< Ce^ (||Vq £ || 2 4 + ||qAq e || 2 ) ds + | jf' ||V • (A^)|| 2 



< (Hql 2 ^ + llVq 5 !!^) (||Vq £ || 2 + ||Aq £ || 2 ) 



ds 



+| J ||V-(A^)fds 
£ (" : + i / ||V-(AV> £ )|| 2 <fe, 



and 



K 9 <e f || V • (Aq e )|| 2 ds + E /"* ||V • (A^ E )|| 2 
io Jo 



K w < f \\V{A9 £ )\\ 2 ds + C f \\A^ £ p\\ 2 ds + C f \\Vi; £ ■ (Vp)\\ 2 ds 
Jo Jo Jo 

+C f \\V 2 p- ip £ \\ 2 ds + C [ || Aq# £ || 2 (is 
Jo Jo 



ft ft 
+C ||V 2 £ -q|| ds + C / ||Vq • (V6 £ )\\ 2 ds 
Jo Jo 



< fwViAO^fds + ^Klo. 
Jo i=i 

Here K\ Q — Kf are estimated as follows: 



Kl + K 2 w <C J (||fl&» + ||Vp||£.) (||VV/|| 2 + ||A^|| 2 )ds 

'||V^ e || 2 + ||A^ £ || 2 ) ds, 
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< C 



2 , Ilv72~||2 

m + Il v P\\m 



(3.21) 



(3.22) 
(3.23) 



(3.24) 



(3.25) 



Kf + Kf <C J {\WW\oo + WWlJ) (\\Apf + ||Aq|| 2 ) ds 

~ c L (" V ^ e "^ 1 + ll WE II^O (ii a pII 2 + ll A qll 2 ) ds ( 3 - 26 ) 

and 

Ki + Kf <C + [| VqHioc) (||Vrf + || A0 e || 2 ) ds 

<c f (livqii 2 ,, + ||v 2 q ||^) (||vr|| 2 + ||Ar|| 2 ) da. (3.27) 

Plugging the estimates for Kg ~ K\q into (|3.2ip . we get 

|| AVI 2 + l|A0 e || 2 + jf (e ||V • (AV e )|| 2 + ||V 3 rf) ds 
<Ce + cJ* (\\Vp\\ 2 H2 + ||Vqfo) (||^ B fo + WWl*) ds. (3.28) 
Combination of (l3~TTi (1349]) and (13~28|) yields 

ll^ll^ + Pin* + J (e l|VV E ||^ + II V0 e |& 2 ) 
<Ce + cJ* (\\Vp\\ 2 H 2 + ||Vq||^ 2 ) (ll^H^ + Uril^) da. (3.29) 

which, together with Theorem 11.11 and Gronwalls inequality gives (|3.ip . This completes the 
proof of Lemma 13.11 □ 

From Lemma [3. 11 we get (|1.15p . Using Sobolev inequality and (|1.15p . we get (|1.16p . This 
completes the proof of Theorem 11.21 
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